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ABSTRACT 

Long term ( "superorbital" ) periods or modulations have been detected in a wide 
variety of both low and high-mass X-ray binaries at X-ray and optical wavelengths. 
A variety of mechanisms have been proposed to account for the variability properties, 
such as precessing and/or warped accretion discs, amongst others. The All Sky Monitor 
on board the Rossi X-ray Timing Explorer provides the most extensive (~15 years) and 
sensitive X-ray archive for studying such behaviour. It is also clear that such variations 
can be intermittent and/or a function of X-ray spectral state. Consequently, we use a 
time-dependent Dynamic Power Spectrum method to examine how these modulations 
vary with time in 25 X-ray binaries for which superorbital periodicities have been 
previously reported. Our aim is to characterize these periodicities in a completely 
systematic way. Some (such as Her X-l and LMC X-4) are remarkably stable, but 
others show a range of properties, from even longer variability time-scales to quite 
chaotic behaviour. 



1 INTRODUCTION 



The Galactic population of X-ray binaries (XRBs) are X-ray 
luminous, interacting binaries in which a donor star (M2) 
transfers material to a neutron star (NS) or black hole (BH) 
compact object (Mx). Low mass X-ray binaries (LMXBs) 
contain donor stars < IMq and typically have orbital pe- 
riods (Port) ranging from hours to days (Liu, van Paradijs 
& van den Heuvel 2007), while high mass X-ray binaries 
(HMXBs) have donor stars > IOMq and P or b ~ several 
days to tens of days (Liu, van Paradijs & van den Heuvel 
2006). 

From the very first decade of satellite X-ray astron- 
omy (the 1970s), long-term monitoring revealed modulations 
that were quasi-periodic on time-scales substantially longer 
than their well-established orbital periods. These included 
results from missions such as Ariel V's All Sky Monitor 
(e.g. Kaluzienski et al. 1976) and Vela 5B (e.g. Priedhorsky 
& Terrell 1983, 1984). Her X-l, with its 35 d on/off cycle 
(which is remarkably stable, at ~ 20 x P or b) exhibits the 
prototypical "superorbital" period (Petterson 1977). 



1.1 Long-term variations 

Such superorbital variations (P sup ) are seen to occur on 
time-scales of tens to hundreds of days, and are thought 
to be related to the properties of the accretion disc (and 
occasionally linked to the donor). Many of these have been 
determined for a number of sources, using archival data from 
later X-ray satellite missions such as CGRO's BATSE (e.g. 
Robinson et al. 1997), RXTE's All Sky Monitor (ASM; e.g. 
Wen et al. 2006) and Swift's Burst Alert Telescope (e.g. 
Farrell, Barret & Skinner 2009). Most recently, superorbital 



behaviour and the various mechanisms proposed to account 
for it, have been summarised in Charles et al. (2008, 2010). 

Radiation-induced warping (Ogilvie & Dubus 2001, 
hereafter OD01) and/or precession (Whitehurst & King 
1991) of the accretion disc may lead to periodic/quasi- 
periodic superorbital modulation of the X-ray flux, e.g. Cyg 
X-2, Her X-l, LMC X-4 & SMC X-l (Clarkson et al. 2003). 
However, superorbital variations may also occur as the result 
of other mechanisms. Here we briefly introduce the physical 
processes that have been proposed. 



1.1.1 Disc Precession 

Precession of accretion discs may occur if q <0. 25-0. 33, due 
to resonances excited in the discs resulting from tidal in- 
teraction with the donor (Whitehurst & King 1991). This 
occurs in the SU UMa sub-class of Cataclysmic Variables 
(CVs; Warner 1995) and the soft X-ray Transients (SXTs; 
O'Donoghue & Charles 1996), containing white dwarfs 
(WD) and BH compact objects respectively. It is very likely 
also to occur in ultra-compact X-ray binaries (UCXBs; e.g. 
X1916-053) where the donor itself is degenerate (Nelemans 
2004). 



1.1.2 Radiation-induced Warping /Tilting 

This is the original superorbital mechanism that was pro- 
posed to account for the 35d cycle in Her X-l by Petter- 
son (1977). Subsequently developed by Wijers and Pringlc 
(1999), the most detailed analysis to date is that of OD01. 
They showed how accretion discs can develop radiation- 
induced warps that are functions of their binary separation 
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Figure 1. OD01 stability analysis for radiation-driven warping 
in X-ray binaries, using a = 0.3 and e = 0.1 (adapted from 
Fig. 7 in OD01, including only sources with known P S up), where 
r;, has units [GMj/c 2 ] and q is dimensionless. The upper two 
lines (solid & dashed) correspond to the first two bending modes 
for r (1£ j ( j = r c (mode starting higher than mode 1). Sustained 
stable warps are only possible close to the solid portions of those 
curves, while warps become increasingly variable/unstable as the 
dashed curves are approached and crossed. The instability zone 
lies between these solid and dashed curves. The bottom solid line 
corresponds to mode 1 for = r , below which discs are un- 
likely to experience radiation-driven warping. The region between 
the bottom and top solid lines are considered an intermediate in- 
stability zone where a disc may display warping cycles as it al- 
ternates between being warped and being flat. Squares indicate 
LMXBs and circles HMXBs. The grey vertical dashed line indi- 
cates the q ~0.33 boundary (Whitehurst & King 1991), to the 
left of which sources are susceptible to disc precession resulting 
from tidal interactions with the donor. 

(rt) and mass ratio (q — K'h/Mx). The stability of such 
warps is defined in this analysis by the range of values for 
the mass input radius (r a dd), with boundary values set by 
the circularization radius (r c ) and the outer disc radius (r ). 
(See FigfT]for further details.) 

OD01 ignored tidal torque and assumed a viscosity pa- 
rameter (q = 0.3) and accretion efficiency (e = 0.1) that 
were considered reasonable for persistent X-ray binaries in 
general, ft is expected that chaotic (e.g. Cyg X-2) or un- 
stable (e.g. SMC X-l) warping may result in quasi-periodic 
superorbital variations, while precession of stable/sustained 
warps (e.g. Her X-l & LMC X-4) should result in steady, 
almost strictly periodic superorbital modulations. 



1.1.3 Magnetic Warping 

Pfeiffer & Lai (2004) carried out numerical simulations of 
viscous accretion discs driven by the magnetic torque in the 
inner disc, resulting from the interaction between the in- 
duced electric currents in the disc and the accreting object's 
magnetic field (£^). They show that an accretion disc sur- 
rounding a magnetic star (T Tauri star, WD or NS) can 
develop warps in the inner disc that steadily precess. The 
retrograde precession as observed in Her X-l, is successfully 
reproduced by their simulations. They suggest that changes 
in M can cause discs to switch between warped and flat. 

1.1.4 Wind-driven Tilting 

Schandl & Meyer (1994) proposed an alternative mechanism 
for the tilted disc observed for Her X-l, involving the coronal 
wind torque in combination with viscous and tidal torques, 
wherein the coronal wind exerts a repulsive force on the disc 
(but see also the analysis and discussion concerning the coro- 
nal model in Scott, Leahy & Wilson (2000)). Quillen (2001) 
proposed that a supersonic wind passing over the accretion 
disc may cause surface instabilities. Moreover, Montgomery 
& Martin (2010) suggest it may also cause lift due to the 
differing speeds above and below the disc, causing a pitch of 
the disc around its line of nodes, thereby tilting it. 

1.1.5 X-ray State Changes 

Long-term modulation of the mass transfer rate (M) 
through their discs has been associated with dramatic "state 
changes" in the X-ray spectra of many sources (of which 
Cyg X-l and GX339-4 have provided frequent demonstra- 
tions). This is particularly true for the SXTs (McClintock & 
Remillard 2006 and references therein) and luminous galac- 
tic bulge sources (van der Klis 2006 and references therein). 
In such cases, the mass transfer rate is intermediate between 
the stable high and low spectral states, and so occasional 
transitions will occur. Super-orbital periods detected in such 
systems may therefore indicate the time-scales associated 
with such transitions (see e.g. King et al. 1997). 

Quasi-persistent XRBs (e.g. KS 1731-260 & X1636-536) 
can display prolonged high (on) and low (off) intervals, char- 
acterized by markedly different flux levels during those times 
(Charles et al. 2008). During the decline from the high to the 
low state, X1636-536 displayed an anti-correlation between 
the hard and soft X-rays as well as superorbital variations 
consistent with modulations of M comparable with the vis- 
cous time-scale in the outer disc (Shih et al. 2005, 2011). 

1.1.6 Precessing Relativistic Jets 

The archetypal galactic micro-quasar is the well-known 
SS433 (Hjellming & Johnston 1981; Margon 1984), which 
displays the remarkable "moving lines" in its optical spec- 
trum with radial velocities that are truly relativistic. These 
are now well-established as arising from precessing relativis- 
tic jets, a process which imparts substantial stability to the 
162 d superorbital period. The jets and their corkscrew mo- 
tion have been directly (spatially) observed at radio wave- 
lengths (Blundell, Bowler & Schmidtobreick 2007), although 
the detailed properties of the 13 d orbital period XRB at 
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the heart of SS433 remain poorly determined (e.g. Barnes 
et al. 2006) . The driving mechanism for the jet precession in 
SS433 is still a matter of controversy, but irradiation driven 
outflows are likely to play a significant role, as discussed in 
e.g. Begelman et al (2006). 

1.1.7 Third Body 

The most luminous LMXB in a galactic globular cluster is 
X1820-303 (at the core of NGC 6624). With close to the 
shortest orbital period (~12 mins), it should not be suscepti- 
ble to radiation-induced disc warping as suggested by ODOf . 
But its very low q (the donor is a very low mass degenerate 
helium dwarf), does make it a prime candidate for tidal disc 
precession effects. Nevertheless, the extremely stable super- 
orbital period which has been observed in X1820-303 is con- 
sidered a likely consequence of the effect of a third body on 
the M in the system (Chou & Grindlay, 2001). Such triple 
systems involving an XRB have long been mooted, especially 
within the extremely dense stellar regions of globular cluster 
cores, but this system remains the strongest candidate. 

1.1.8 Be X-ray Binaries 

When a NS in a highly eccentric orbit passes through, or 
interacts with, the varying equatorial disc surrounding a 
rapidly rotating Be-star, it naturally gives rise to a modula- 
tion at the orbital frequency (see e.g. Leahy & Kostka 2008 
and references therein, and Hayasaki & Okazaki 2006) . How- 
ever, there are longer term variations in which the equatorial 
disc expands and contracts which will give rise to a modula- 
tion of the accreted fraction on this superorbital time-scale 
(Charles et al. 2008). These are the Be X-ray binaries (BeX), 
which display a wide variety of variations at X-ray & opti- 
cal wavelengths, extending to time-scales much longer than 
being considered here (Rajoelimanana, Charles & Udalski 
2011). 

1.2 Time-dependent Period Analysis 

Periods are typically identified by applying a periodic analy- 
sis to the entire available dataset to produce a periodogram. 
While this approach allows one to distinguish between peri- 
odic (sharp peaks) and quasi-periodic (broadened or multi- 
ple peaks) behaviour, it is essentially limited to the identifi- 
cation (not characterization) of periodic signals (Wen et al. 
2006). 

However, from our brief introduction to the variety of 
physical mechanisms that can induce superorbital varia- 
tions, it will be clear that some of these can be intermit- 
tent, or vary in time-scale so much as to be rendered un- 
detectable in a summed periodogram. Consequently, such 
behaviour should be analysed in a way that displays the 
"instantaneous" modulations present in the data, and then 
follow this as a function of time. 

A time-dependent periodic analysis of two HMXBs and 
two LMXBs in the form of Dynamic Power Spectra (DPS; 
e.g. Clarkson et al. 2003), clearly illustrated the need for 
a more comprehensive approach in the study of superor- 
bital periods. Therein the P 3up in SMC X-l was seen to 
evolve dramatically and systematically, while those in Her 



X-l and LMC X-4 remained essentially steady in period, 
although the modulation in Her X-l would occasionally dis- 
appear completely e.g. Still et al. (2004). In contrast, the 
Psup observed in Cyg X-2 appeared to be highly variable 
and unstable. 

The DPS technique is well established and the X-ray 
burst community uses a variant of it routinely (e.g. Fig. 19 in 
Watts et al. (2005)). However, since it requires significantly 
more computations than the typical period analysis, its ap- 
plication was limited in the past. Recent advances in desktop 
technology allow this technique to be employed more readily 
than before. 

Another approach to unravel the time-dependence of 
these long-term periods was recently developed by Hu et 
al. 2011, and this is the Hilbert-Huang Transform (HHT). 
Hu et al applied the HHT to the same SMC X-l dataset as 
used here, and find that the HHT is consistent with results 
obtained from the DPS. They suggest that it may be able 
to show more detail in the frequency and time domains. 

Boyd & Smale (2004) applied time-frequency decompo- 
sition techniques to the long-term RXTE ASM lightcurves 
of three XRBs, separating them into random and periodic 
components. They suggest that all the P sup modulations 
in Cyg X-2 are integer multiples of its P or b, although they 
concede that competing warping modes (suggested in OD01) 
could be considered as one of several mechanisms which may 
possibly be responsible for its long-term variations. 



1.3 Target Details 

With an RXTE observational baseline that has doubled, we 
revisit the 4 sources of Clarkson et al. (2003). In addition, 
we apply the same analysis to the other 14 sources showing 
Psup behaviour and summarised in Charles et al. (2008). A 
further 7 sources, not contained therein or published since, 
were also included in our target list. 

Here we present a time-resolved periodic analysis of all 
25 sources with previously published or claimed P sup <1 
yr, many of which have been inferred to be associated with 
warped and/or precessing accretion discs. Consequently, we 
endeavour to provide a comprehensive overview of the long- 
term behaviour of the majority of superorbital modulations 
observed in X-ray binaries, through this systematic DPS ap- 
proach. In this way we hope that these results can provide 
a fresh context within which superorbital periods in X-ray 
binaries can be investigated. 

The 9 HMXBs and 16 LMXBs considered here are listed 
in Tables [TJ & [2] respectively, together with their previously 
published P sup , P or b and q — jj 2 - (values of which are 
taken from OD01 unless otherwise stated). Their positions 
(if known) are shown on an adapted OD01 plot (Fig[T]) of 
instability zones with respect to ri, (in units of [GMx /c 2 ]) 
and q. 

2 OBSERVATIONS AND DATA ANALYSIS 
2.1 RXTE/ ASM 

The Massachusetts Institute of Technology (MIT) has oper- 
ated an All Sky Monitor (ASM) on board the Rossi X-ray 
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Table 1. HMXBs with known supraorbital periods 



Source 


P S u P [d] 


P r6 [d] 


3 


Cen X-3 


140 1°) 


2.09 W 


17.0 


Cyg X-l 


142 326 M 


5.6 M 


1.7 




150 & 290 l ia ' 






LMC X-3 


99 [ C2 1, 100-500 W 


1.70 


0.5 


LMC X-4 


30 W 


1.41 


10.6 


SMC X-l 


50-70 W 


3.89 l s l 


11.0 


SS433 


162 W 


13.10 [°1 


[1.0] 


X01 14+650 


31 W 


11.6 W 




X1907+097 


42 [ p l 


8.38 W 


12.0 


XTE J1716-389 


99 W 







I 01 ' Cowley ct al. (1983), [ C2 i Cowley et al. (1991), 
' c l Crampton, Cowley & Hutchings (1980), 
OD01 is the only reference found for it, 
l La i Lachowicz ct al. (2006), 

' L l Liu, van Paradijs & van den Heuvel (2006) and ref. therein, 
l p l Pricdhorsky & Terrell (1984), 

Rico (2008), 
' s l Schreier et al. (1972), 
W Wen et al. (2006) and ref. therein 



Table 2. LMXBs with known superorbital periods 


Source 


Psup [d] 


^orb 


[d] 


9 


Cyg X-2 


50-80 1°] , 60-90 l w ) 


9.844 


[L] 


0.34 


EXO 0748-676 


181 M 


0.158 


[L] 




GRS 1747-312 


147 [W] 


0.515 


m 




GX 339-4 


190-250 ' c l 


1.755 


m 


[0.7] 


GX 354-0 


63 or 72 W 








Her X-l 


33-37 [ ie l 


1.700 


w 


1.56 


ICR J 17098-3628 


163 W 








KS 1731-260 


38 








LMC X-2 


8 Id 


0.34 


m 




MS 1603.6+2600 


5 1°] 


0.077 






Sco X-l 


2.6 * 


0.788 


w 


[0.7] 


X1636-536 


46 Id 


0.158 




[0.6] 


X1730-333 


217 W 








X1820-303 


172 [^1 


0.008 


w 


[0.1] 


X1916-053 


5 [C] & 199 l p l 


0.035 




[0.1] 


X1957+115 


117 [°1 


0.390 


w 


[0.7] 



[ c l Charles et al. (2008) and ref. therein, 
W Kotze, Charles & Crause (2009), 

I L 1 Liu, van Paradijs & van den Heuvel (2007) and ref. therein, 

I Le l Leahy & Igna (2010), 

[ p l Priedhorsky & Terrell (1984), 

t w l Wen et al. (2006) and ref. therein, 

* Quoted as ~62 d (OD01), not ~62 hrs (Kudryavtsev et al. 1989) 

Timing Explorer (RXTE) since early 1996. The ASM ob- 
serves the X-ray sky by using 3 rotating Scanning Shadow 
Cameras (SSC) to scan ~80% of the sky during each ~90 
min orbit. Data are reduced and compiled weekly by the 
ASM team and made publicly available as dwell-by-dwell 
or one-day-averages in four energy bands, which includes 

1 |http://xte.mit.edu/|A SMlc.html 



a Sum-band (1.5-12 keV). One-day-average data are the 
dwell-by-dwell data binned into 1-day bins. Archival ASM 
datasets contain the lightcurves of all X-ray sources in the 
RXTE catalogue. For full details see Levine et al. (1996) and 
http:/ /xte. mit.edu/. 

2.2 DPS Variability Analysis 

The ASM one-day-average datasets used in this paper span 
from 20 February 1996 to 12 February 2011 (MJD 55608) 
and therefore provide 15 year light-curves of unparalleled 
quality and sensitivity. While period analysis of the entire 
dataset can provide the maximum sensitivity to low-level 
modulations (if they are steady), it has serious limitations 
when dealing with unstable or evolving periodic signals. Sig- 
nals that are intermittent may be completely damped out, 
and evolving signals may cause multiple peaks over a much 
larger frequency range than is expected for so-called quasi- 
periodic oscillations (see e.g. van der Klis 2006). These ef- 
fects get worse with longer observational baselines. Since 
superorbital periods relating to accretion disc behaviour are 
expected to be quasi-periodic, or could be unstable and may 
drift systematically, our interest rather lies in characterising 
the behaviour of those periods over time. 

2. 2. 1 Windowing 

The Dynamic Power Spectrum (DPS) method employed 
by Clarkson et al. (2003) has the advantage of illustrat- 
ing very clearly whether periodic signals are sustained, and 
whether they are stable or drifting in period (although the 
method is also sensitive to phase drifts; Clarkson et al. 2003). 
The method requires the datasets to be split into windows 
that are of sufficient length to allow detection of the maxi- 
mum period considered. A Lomb-Scargle (L-S) periodogram 
(Lomb 1976 & Scargle 1982, 1989) was produced for every 
such dataset and the results of all the periodograms for a 
source were plotted together in a density map, using the L- 
S power for each frequency plotted at every window's mid- 
point along the time-axis. The frequency domain covered 
periods with a minimum of 2 d to a maximum comparable 
to the window size. 

Larger window sizes cause the smearing out of variabil- 
ity as it averages out small variations, but it has the ad- 
vantage of enhancing sustained periodic signals. However, 
variable or intermittent periodic signals may be completely 
washed out if windows are too large. Window sizes were ini- 
tially chosen to be 400 d for all sources, since it provides 
coverage of ~5-10 periodic cycles in the range ~40-80 d and 
it has the advantage of being directly comparable to the ear- 
lier results already published for SMC X-l, LMC X-4, Her 
X-l & Cyg X-2 (Clarkson et al. 2003). It also allowed us to 
identify longer potential periods, which could then be inves- 
tigated in more detail by choosing more appropriate window 
sizes to sample them properly. For periods >80 d, windows 
of ~5 times the periods were used for the results presented 
here, and for periods <10 d, we used 100 d windows. 

A sliding window approach was also employed, whereby 
consecutive datasets overlap such that they move in 50 d 
steps in the time domain. This approach provides adequate 
temporal resolution, whilst smoothing out the noise from 
otherwise independent windows that will occur. 
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2.2.2 Dynamic Window Function 

To ensure that the features contained in the DPS are not ar- 
tificially induced by the sampling or windowing employed, 
the spectral window function was determined for each 
dataset and a density map was constructed with that infor- 
mation. The resulting Dynamic Window Function (DWF) 
was compared to the DPS for each source. No significant 
features contained in the DPS were repeated in the DWF of 
any source considered here. 

2.2.3 Dwell-by- dwell data 

The DPS analysis was also repeated for 1-day averages that 
were constructed from the dwell-by-dwell data, employing 
aggressive filtering of the data to eliminate possibly dubious 
data points, as suggested on the RXTE ASM website. Those 
results are consistent with what we present here. However, 
the data loss due to the more conservative approach, does 
lead to diminished temporal coverage. We find that the fea- 
tures in the DPS are in general clearer using the ASM 1-day 
averages, which allows better coverage, and are therefore the 
results presented here. 

2.2.4 Noise Levels 

White noise (frequency independent) levels were determined 
for each complete 15 year one-day-average dataset by gen- 
erating 10 4 random datasets by Monte Carlo simulation, as- 
suming a white noise distribution and using the time values, 
mean and standard deviation of the dataset. Lomb-Scargle 
periodograms for each random dataset yield periods associ- 
ated with the maximum power in each, and the subsequent 
distribution of powers amount to a cumulative probability 
distribution function, providing the false alarm probability 
(result of white noise) associated with a particular power. 
The white noise level for a 99.99% confidence level, is the 
power associated with a false alarm probability of 0.01%. 
The noise level for the dataset as a whole can be used as 
an estimate for its windows and is not expected to differ 
significantly between sources. 

Red noise (frequency dependent) levels for the RXTE 
ASM data, determined by Farrell, Barret & Skinner (2009) 
for 5 of these sources, indicated that the red noise levels do 
not always exceed the white noise at the low frequencies. For 
Sco X-l and X1916-053 the white noise levels were ~doublc 
those of the red noise, while Cyg X-2, X1636-536 & X1820- 
303 had red noise levels that rose to >3 times the white noise 
levels at low frequencies (0.001-0.01). Apart from differing 
significantly for each source, these modelled red noise levels 
are extremely dependent on the time domain covered by the 
dataset, as well as the weighting scheme applied to it, if any. 

However, the aim of this paper is to characterize the 
behaviour of previously reported periodic signals over time, 
not claim the discovery of new periodic behaviour. The latter 
would normally require rigorous testing against noise, while 
such requirements should be irrelevant in case of the for- 
mer. While noise levels may serve as an accepted statistical 
measure for the significance of prominent periodic signals, 
the true measure of their significance lies in their persis- 
tence over time, which a time-dependent periodic analysis 
provides automatically. 



2.3 Results 

To allow easier interpretation of features in the power den- 
sity maps (the 2-D representation of 3-D data) and how they 
translate to features contained in the lightcurves, the DPS 
(from the Sum-band one-day-average data) are plotted, to- 
gether with the lightcurve of each source. The light-curves 
shown were re-binned into 10-day bins and plotted above 
the DPS, over the same temporal range (horizontal/x-axis). 

Furthermore, the L-S periodogram for the entire 15 year 
dataset was plotted to the left of the DPS for each source 
and over the same frequency range (vertical/y-axis). This 
allows immediate comparison of features in the DPS to the 
overall periodogram. Frequency labels were replaced by the 
equivalent period labels, to facilitate easier comparison with 
previously published P sup - 

Known P sup contained in Tables [T]& [2] are mentioned in 
the caption of each source and are also clearly indicated on 
each plot with ticks along the period/frequency-axis (red in 
the online version of this paper). The density scale-bar in- 
cluded in each plot indicates the power (line-of-sight/z-axis) 
associated with the periodogram, with stronger detections 
being darker. 

DPS and L-S plots were produced for frequencies ap- 
propriate to their reported superorbital periods. Including 
a wider range of frequencies around those periods allows a 
more complete picture of the temporal properties of each 
source. In an effort to allow comparability between sources, 
the size of the frequency range plotted was kept relatively 
constant at ~0. 02-0. 03 for all the plots, except that sources 
with Psup <10 d were plotted over a frequency range of 0.1. 
We only include the plot for the P sup ~199 d in X1916-053 
(Fig |25|l . since we find no evidence for the P aup ~5 d. 

The estimated white noise levels for all sources are, not 
surprisingly, quite similar, being in the range from 15.6-16.8; 
they are plotted as vertical lines (green in the online version 
of this paper) on the L-S plots of the sources as a guide. 
Prominent features in the L-S periodogram which are well 
above these estimated white noise levels, can normally be 
considered significant, the situation for the majority of the 
sources considered here. Only in Figs. [9] [19] & [20] were no 
significant features detected in the L-S or DPS at any time. 

Our aim here is to present the results in a way that pro- 
vides the reader with sufficient information to reach their 
own conclusion regarding the variability properties of a 
source at a single glance. All relevant additional informa- 
tion, not contained in the plots themselves, was therefore 
included in the caption of each source. The figures for the 
sources are presented in the same order as their listing in Ta- 
bles [T]& [2] with the HMXBs presented before the LMXBs. 

While the results presented here certainly include most 
of the X-ray binaries with published superorbital periods, 
we do not claim to have included all sources for which su- 
perorbital periods have been published, nor do we claim to 
have included all X-ray binaries that have ever displayed su- 
perorbital behaviour. However, those included here do have 
superorbital periods that have been proposed to be associ- 
ated with warped and/or precessing accretion discs at some 
stage. Their inclusion is intended as an opportunity to study 
such features using the DPS approach, thereby indicating 
whether such modulations were steady, evolving, persistent 
or intermittent. 
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Figure 2. Cen X-3 : HMXB with P sup ~140 d (in OD01 only), 
indicated by marker on y-axis. Top panel: RXTE ASM lightcurve, 
Left panel: L-S over entire ~15 yr dataset, Main panel: DPS with 
time vs period (x-y), Right: scale bar indicating DPS power (z). 
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Figure 3. Cyg X-l : HMXB with P sup ~142 d (in OD01 only), 
~150 & 290 d (Lachowicz et al. 2006) and ~326 d (Rico 2008). 
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Figure 4. LMC X-3 : HMXB with P sup ~99 d (Cowley et al. 
1991). 
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Figure 5. LMC X-4 : HMXB with P sup ~30 d (Wen et al. 2006). 
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Figure 6. SMC X-l : HMXB with P sup ~50-70 d (Wen et al. 
2006). 
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Figure 7. SS433 : HMXB with P sup -162 d (Wen et al. 2006). 
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Figure 8. X0114+650 : HMXB with P sup ~31 d (Wen et al. 
2006). 
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Figure 11. Cyg X-2 : LMXB with P sup ~50-80 d (Charles 
2008) or 60-90 d (Wen et al. 2006). 
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Figure 9. X1907+097 : HMXB with P 3up ~42 d (Priedhorsky 
& Terrell 1984). All features are below white noise. 
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Figure 12. EXO 0748-676 : LMXB with P sup 
higher state (Kotze, Charles & Crause 2009). 
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Figure 10. XTE J1716-389 : HMXB with P 3up ~99 d (Wen et 
al. 2006). 
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Figure 13. GRS 1747-312 : LMXB with P sup ~147 d (Wen et 
al. 2006). 
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Figure 14. GX 339-4 : LMXB with P sup ~190-250 d (Charles 
et al. 2008). 
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Figure 17. IGR J17098-3628 : LMXB with P s , 



^163 d during 



part of the high state (Kotze, Charles & Crausc 2009). 
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Figure 15. GX 354-0 : LMXB with P 3up ~63 or 72 d (Liu, 
Paradijs & van den Heuvel 2007). 
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Figure 18. KS 1731-260 : LMXB with P aup ~38 d (Charles et 
al. 2008). 
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Figure 16. Her X-l : LMXB with P sup ~33-37 d, with average 
~35 d (Leahy & Igna 2010). 
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Figure 19. LMC X-2 : LMXB with P sup ~8 d (Charles et al. 
2008). All features are below white noise. 
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Figure 20. MS 1603.6+2600 : LMXB with P sup ~5 d (Charles 
et al. 2008). All features are below white noise. 
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23. X1730-333 : LMXB with P sup ~217 d (Charles et al. 
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Figure 21. Sco X-l : LMXB with P sup ~2.6 d (Kudryavtsev et 
al. 1989). 
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Figure 24. X1820-303 : LMXB with P sup ~172 d (Wen et al. 
2006). 
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Figure 22. X1636-536 : LMXB with P sup ~46 d after transition 
to lower state at MJD~52000 (Charles et al. 2008). 
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Figure 25. X1916-053 : LMXB with P sup ~199 d (Wen et al. 
2006). 
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Figure 26. X1957+115 : LMXB with P sup ~117 d (Charles et 
al. 2008). 



& Wilms (2007), now appears to be repeated approximately 
10 years later (around MJD~54000). There also appear to 
be a decrease in maximum and increase in minimum 10- 
day binned flux coincident with the temporal location of the 
minimum superorbital periods in those two features in the 
DPS curve. The behaviour after these two largest features 
does not repeat, but appears to be similar in terms of the 
decreased period range covered, and contains features that 
seem to recur on a time-scale of ~1000 d. 

SMC X-l lies in the instability zone on FigfTJ How- 
ever, it is situated closer to the dashed curve and therefore 
a radiation-induced warp in its accretion disc is expected to 
be variable, though not yet unstable. The steady evolution of 
its superorbital period may subsequently be the result of the 
precession of a steadily warped disc in which the warp itself 
varies, due to the competition of warping modes (Clarkson 
et al. 2003). Full investigation and further interpretation of 
these results warrants a far more detailed discussion of the 
source itself that goes beyond the scope of this paper. 



3 DISCUSSION 

It is only fitting that we begin our discussion with the up- 
to-date results for the 4 sources included in Clarkson et 
al. (2003). We will then discuss the remainder of those in 
Charles et al. (2008) plus 7 additional sources. For each 
source we explicitly discuss its stability to radiation-driven 
warping under the OD01 framework. We remind the reader 
that we are not claiming that radiation-driven warping must 
explain superorbital periodicities for all sources. 

Determining the behaviour of the previously reported 
superorbital periods required a time-dependent method, as 
suggested in OD01. They predicted that systems able to 
produce steady superorbital periods resulting from stable, 
steadily precessing, warped discs (as in Her X-l), should be 
rare, and variable/chaotic warping would result in a variety 
of quasi-periodic or aperiodic behaviour in the majority of 
sources. Our results here show that this prediction is largely 
borne out. While Fig[T] provides us with the OD01 stability 
predictions for individual sources, some were discussed in 
more detail in OD01 and where relevant we refer to those 
comments, since they pertain to deviations from the predic- 
tions implied in the figure. 

3.1 The Clarkson et al. (2003) sources 

Our DPS results for Her X-l, SMC X-l, LMC X-4 and Cyg 
X-2 are consistent with those presented in Clarkson et al. 
(2003) and those included in the Charles et al. (2008) re- 
view. However, the observational baseline used here has ap- 
proximately doubled. These four sources therefore represent 
good examples with which the long-term behaviour of the 
superorbital periods of other sources can be compared, par- 
ticularly in terms of OD01 disc instability. 

3.1.1 SMC X-l (FigUty 

Still by far the most intriguing of all our results, the DPS of 
this HMXB clearly shows the smooth and steady evolution 
of its Psup over the range ~42-70 d. The behaviour early in 
the light-curve (the feature surrounding MJD~50500), also 
included in Clarkson et al. (2003) and Trowbridge, Nowak 



3.1.2 LMC X-4 (Fig\3j) 

Perhaps the steadiest of the stable superorbital periods in 
all HMXBs, a substantial modulation is strongly detected 
with a Psup of 30 d. The detection level appears to drop 
temporarily at MJD~51000, MJD~53000 and finally almost 
disappears completely around MJD~55500. Since the source 
is clearly detected throughout our observational baseline, 
the disappearance of the periodic signal toward the end is 
not due to poor coverage, but more likely an effect intrinsic 
to the source itself. 

LMC X-4 lies in the instability zone on Fig[TJ closer to 
the solid curve. The radiation-induced warps in its accretion 
disc are therefore expected to be stable and consequently 
able to produce a steady modulation. 

3.1.3 Her X-l (Fia YTty 

This prototypical LMXB has continued to show a very 
strong detection of a relatively stable P sup ~35 d through- 
out the ASM coverage (and whenever it has been moni- 
tored over the last 40 years). The interruptions occurring 
around MJD~51250-51750 and MJD~53000, coincide with 
the anomalous low states (ALS), of which the first was in- 
cluded and discussed as such in Still et al. (2001). During 
these states, the maximum flux drops dramatically and be- 
comes comparable to the minimum flux, which also appears 
to be decreasing slightly. Multi-wavelength data suggest that 
the ALS is not due to M being cut off, but is rather caused 
by the shadowing effect of a more extreme tilt in the disc, 
which blocks the X-ray emission from the inner regions from 
the observer (Vrtilek et al. 2001). Leahy & Dupuis (2010) 
suggest that a larger disc warp may account for the larger 
reduction in extreme UV flux (using Extreme Ultraviolet 
Explorer data) in comparison to that in X-rays (using si- 
multaneous RXTE ASM data) during the ALS. 

Her X-l lies in the instability zone on Fig[TJ approx- 
imately halfway between the solid and dashed curves. The 
superorbital period is therefore expected to be relatively sta- 
ble, as a result of the precession of a radiatively warped ac- 
cretion disc. However, it may show some evolution due to 
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small-scale variation in the structure of the warp itself (see 
Leahy & Igna 2010). 

3.1.4 Cyg X-2 (Fig\T7j) 

Still showing chaotic and unstable behaviour, the range of 
peaks from ~40-90 d in this LMXB are strongly detected 
from time to time throughout the RXTE lifetime. An ini- 
tially weaker (but significant) detection in the range ~40-50 
d becomes a much stronger detection toward the end of the 
light-curve. 

Cyg X-2 lies well above the dashed curve on Fig[TJ and 
is expected to experience chaotic warping of its accretion 
disc. As such, it is not expected to have a stable or steadily 
evolving superorbital period, but rather intermittent and 
variable periods produced by precessing unstable warps. The 
DPS is consistent with intermittent variations. 

3.2 Other HMXBs 

3.2.1 Cen X-3 (Fig® 

There is a detection of the P sup ~140 d around MJD~51500 
for this HMXB, but it seems to form part of an evolving fea- 
ture in the DPS, that had its strongest detection at ~200 
d at MJD~51000. A variety of shorter term signals (proba- 
bly linked to the time-scale of the "on" states) are strongly 
detected up to MJD~53500, after which the maximum flux 
drops dramatically. 

Cen X-3 is located in the OD01 instability zone in Fig[TJ 
between LMC X-4 and SMC X-l. Consequently, Cen X-3 is 
expected to develop a warped precessing accretion disc in 
which the warp itself might be variable (Priedhorsky & Ter- 
rell 1983; Iping & Petterson 1990). The behaviour of the 
superorbital period is therefore also expected to be some- 
where between stable (like LMC X-4) and steadily evolving 
(like SMC X-l). However, our DPS shows the evolution to 
be erratic, but not chaotic (like Cyg X-2). 

3.2.2 Cyg X-l (Fig® 

In this HMXB we detect the very long term P sup of ~150 
d and ~290 d, which were detected in the Ariel 5 and Vela 
5B datasets respectively (Lachowicz et al. 2006), and the 
subsequently detected ~326 d period in the Swift BAT and 
RXTE ASM data (Rico 2008). In fact, apart from those we 
detect several other features, many of which appear to be 
stable and are not necessarily harmonics. 

Cyg X-l is located on the OD01 curve associated with 
r a dd = r c on FigfT] where we expect stable warps to de- 
velop resulting in a steady superorbital period, like that in 
LMC X-4. Zdziarski, Pooley & Skinner (2011) reported the 
doubling of the P S u P and the constant shape of the intrinsic 
spectrum, concluding that those modulations are likely pro- 
duced by a precessing accretion disc and not changes in M. 
They also confirm the presence of a bulge on the accretion 
disc edge, as previously proposed. 

3.2.3 LMC X-3 (Fig[§ 

Features in the DPS of this HMXB are closely linked to 
details visible in the lightcurve, since the range of superor- 



bital periods and the clearly visible lightcurve features are 
all ~ 100-500 d. The periodic signals are therefore rather in- 
dicative of time-scales than periodicity (OD01). 

There is an initial evolution from shorter to longer term 
features, then a change of direction around MJD~51500 
which continues past MJD~53000 and then collapses. How- 
ever, a new long term cycle then starts, evolves more rapidly 
to a shorter term feature during MJD~53000-54500 and also 
appears to collapse, making way for another long term cycle 
starting at MJD~54700. We see no evidence for a steady 
Psup ~99 d as reported in the past (Cowley et al. 1991), but 
rather several steadily evolving long-term features covering 
the range ~ 100-500 d, which is the range reported in Wen 
at al. (2006). 

LMC X-3 is located on the OD01 curve using r a dd = r o 
in Fig[TJ below which systems are not expected to develop 
warped precessing accretion discs. However, OD01 suggest 
that the long-term variations in this source are due to 
changes in the mass accretion rate from the donor. 

3.2.4 SS433 (Fig^ 

The strongest detection for this HMXB is P aup ~162 
d, which is the well-known relativistic jet precession pe- 
riod (Margon et al. 1984). This disappears shortly after 
MJD~51000 and reappears at MJD~52000 at a weaker de- 
tection strength, remaining until MJD~55000. The weaker 
detections at the same frequency over the entire baseline and 
the singular strong detection in the overall L-S, suggest that 
a stable X-ray period ~162 d is definitely present. The de- 
tection thereof may just be hampered from time to time by 
intrinsic fluctuations. A strong detection of a period ~400 
d at MJD~51000-52000, when the detection of the ~162 d 
period effectively disappears, makes a weaker reappearance 
at MJD~53000-54000. During that time it is mirrored (on 
the opposite side of the ~162 d period) by an equally weakly 
detected ~100 d period, which appears to be unstable but 
persists to MJD~55000. Furthermore, upon its reappear- 
ance at MJD~51000-52000, the ~400 d period appears to 
have evolved to an even longer period. Since the ~400 d pe- 
riod is detected more strongly when the ~162 d period is 
weakly detected, this might suggest that the presence of the 
longer period may have something to do with the weaker 
detections of the seemingly stable ~162 d period. 

SS433 is located in the instability zone on FigfT] next 
to Her X-l. A stable superorbital period is expected and 
is very clearly present in its optical and radio behaviour. 
The relativistic jets as inferred from their optical emission 
components, have been stable on the 162 d period for > 30 
yrs (Blundell & Bowler 2005). We note that SS433 is in a 
very different evolutionary state at present compared to any 
other XRB considered here, with its high q value leading to 
extreme super-Eddington accretion onto the BH (Begelman, 
King & Pringle 2006). 

3.2.5 X01 14+650 (FigWS) 

In this HMXB, the P SU p ~31 d is weakly detected from 
MJD~50500 in the DPS. It is strongly detected in the over- 
all L-S periodogram and from MJD~50700-51400 in the 
DPS. Thereafter it remains relatively stable with alternating 
stronger and weaker detections of this signal. 
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Although it is not as strongly detected as those of LMC 
X-4 or Her X-l, it appears relatively steady and persistent. 
The location of X0114+650 on FigfT] is n °t known, but as 
an HMXB with P or b ~ 11.6 d (Liu, van Paradijs & van den 
Heuvel 2006) it is likely to be above the OD01 curve for 
r a dd = t , and stable superorbital periods associated with 
a stable precessing warped accretion disc are therefore not 
unexpected. Moreover, since X0114+650 is a BeX (Koenigs- 
berger et al. 1983), we could also expect superorbital varia- 
tions due to a modulation in M as the equatorial disc of the 
Be-star expands and contracts (see Rajoelimanana, Charles 
& Udalski 2011). 

3.2.6 X1907+097 (Fig\§ 

We do not find evidence for the claimed (OD01) P S u P ~42 
d (or any other) X-ray superorbital period in this HMXB. 
X1907+097 is located above the dashed curve on FigfTJwhere 
we expect chaotic warping to occur, but OD01 suggested 
that it may lack superorbital variations as a result of its 
eccentric orbit (Bildsten et al. 1997). 

3.2.7 XTE J1716-389 (Fig\W\j 

XTE J1716-389 is a quasi-persistent HMXB. An initial de- 
tection of the P sup ~99 d, becomes progressively weaker and 
is only visible in the DPS until MJD~53000, when there is 
a dramatic drop in the lightcurve to the low state. 

The relatively regular, weaker detections of the ~99 d 
period while the source is persistent, do seem to suggest 
that it exhibits a relatively steady superorbital period. The 
location of XTE J1716-389 on FigfTJis not known, but as an 
HMXB it is expected to be at least above the OD01 curve for 
T a dd = r . As such, a stable superorbital period associated 
with a steadily precessing warp while it is persistent, is also 
not unexpected. 

3.3 Other LMXBs 

3.3.1 EXO 0748-676 (Fig{7j$) 

A Psup of 181 d was detected during the high state in 
this quasi-persistent LMXB and was clearly not stable, but 
shows steady evolution. The DPS shows several significant 
detected signals during the portions of the higher state 
where the previously published superorbital periods were 
noted (Kotze, Charles & Crause 2009) and there is no evi- 
dence for any periodic signals in the low state. 

EXO 0748-676 has an unknown location on FigfT] How- 
ever, it has P or b ~3.8 hours (Liu, van Paradijs & van den 
Heuvel 2007) and q ~0. 11-0. 28 (Munoz-Darias et al. 2009), 
implying M2 <0.42. As a LMXB with P or b <1 day, we would 
not expect it to produce a steadily precessing warped accre- 
tion disc, but with a q <0. 25-0. 33, tidal disc precession is a 
serious possibility. 

3.3.2 GRS 1747-312 (FigUM 

The P aU p ~147 d of this LMXB is strongly detected in the 
L-S periodogram and DPS at MJD~50500 & MJD~52000. 
However, this periodic signal shows clear variation that 



appears to continue through an interval with lower am- 
plitude variations (significantly weaker detections around 
MJD~52500-53500), to re-emerge thereafter. 

There is also a significant detection at ~70 d, which 
appears to be related to (but not simply the harmonic of) 
the periodic behaviour around ~147 d. The two periodic 
signals at ~70 and ~147 d both appear to evolve overall 
toward slightly longer periods. 

The location of GRS 1747-312 on Fig[TJ is not known, 
but as a LMXB with P or b ~12.36 hours (Liu, van Paradijs 
& van den Heuvel 2007) it is likely situated below the OD01 
curve for r a dd = r , below which sources are not expected to 
develop warped precessing discs or display the steady P 3up 
associated therewith. 

3.3.3 GX 339-4 (Fig\iJ§ 

The P SU p ~ 190-250 d in this transient-like BH-LMXB are 
detected in the lightcurve after removing the outbursts (the 
plot included herein) as well as when the outbursts are in- 
cluded. GX 339-4 is located below the curve for r a£ j«j = r on 
FigfTJand is not expected to develop a warped precessing disc 
and variations are likely due to modulations in M (OD01). 
However, the q ~0.7 used in OD01 was a rough estimate 
and it has since been determined to be q <0.08 (Hynes et 
al. 2003), making tidal disc precession likely. Furthermore, 
Hynes et al. (2003) find P orb ~1.755 d and M x ~5M S . 
Recalculation using the OD01 method gives ~1 in com- 
parison with ~0.6 in OD01, moving GX 339-4 into the in- 
termediate instability zone. 

3.3.4 GX 354-0 (Fig\T$j) 

The DPS shows a Psup ~63 & ~72 d are both initially de- 
tected, and rapidly replaced by a single period between those 
values. The strongest detections are towards the end of the 
dataset where a ~50 d period evolves toward ~70 d. The 
source also clearly displays very long-term variability on a 
~decade long time-scale (Ozdemir 2010, Kotze & Charles 
2010), which are not believed to be associated with the ac- 
cretion disc properties. 

We find no evidence for a stable superorbital period, nor 
for prolonged steady evolution. The behaviour rather resem- 
bles the chaotic warping observed in Cyg X-2. Its location 
on FigfT] is unknown. 

3.3.5 IGR J 17098-3628 (Fig\F§ 

A Psup of 163 d was detected during the higher state in this 
quasi-persistent LMXB. It is clearly not stable and shows 
steady evolution with several significant detected signals 
in the DPS during the portions of the higher state where 
the previously published superorbital periods were noted 
(Kotze, Charles & Crause 2009). There is no evidence for 
periodic signals in the low state. 

IGR J 17098-3628 has an unknown location on FigfT] but 
in Kotze, Charles & Crause (2009) we argued that, based on 
the apparent similarities with EXO 0748-676, IGR J17098- 
3628 likely also has P orb <1 d. As a LMXB with P orb <1 
day, IGR J17098-3628 would not be expected to produce 
steadily precessing warped accretion discs (OD01). 
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3.3.6 KS 1731-260 (Ftg\Wjj 

The Psup ~38 d of this quasi-persistent LMXB is detected 
during the high state (MJD<51000) and during the tran- 
sition from the high state to the low state (MJD~51000- 
52000). However, the strongest detection is a period that 
evolves from ~70 d to ~45 d during that transition. So 
there appears to be two sections to the high state which 
could be considered separately, namely MJD<5f000 and 
MJD~51000-52000. The ~38 d superorbital period appears 
in both of those, suggesting some measure of stability. 

Although KS 173f-260 does not appear to be stable 
like Her X-l when it is persistent, its behaviour is not quite 
comparable to the chaotic warping behaviour observed in 
Cyg X-2. KS 1731-260's location on Fig[T] is unknown, but 
as an LMXB it is likely located below the OD01 curve for 
i~ a dd — f'o and therefore unlikely to produce a warped accre- 
tion disc. 

3.3.7 LMC X-2 & MS 1603.6+2600 (FigUM &Wj) 

The DPS results for these LMXBs show no evidence for 
any significant long-term X-ray periodicities, particularly 
excluding the previously suggested P sup ~8 d and ~5 d for 
LMC X-2 & MS 1603.6+2600 respectively (Charles et al. 
2008). All points in their L-S periodograms are below the 
white noise level and their DPS plots show no periodic fea- 
tures that are persistent in any way. We note that periods 
in the same range are very clearly detected for several other 
sources using similar data. 

The behaviour for LMC X-2 and MS 1603.6+2600 is 
not even reminiscent of Cyg X-2, where there are significant, 
but chaotic signals that are persistent and strongly detected. 
The results for these two sources appear to represent only 
noise. 

Locations for both these sources are unknown on Fig[T] 
MS 1603.6+2600 has P orb ~1.85 hours and LMC X-2 has 
Port ~8.16 hours (Liu, van Paradijs & van den Heuvel 2007) 
and as LMXBs with P or t <1 d, they are not expected to pro- 
duce warped accretion discs (OD01). However, with q ~0.14 
for LMC X-2 (Cornelisse et al. 2007), tidal disc precession 
is likely. MS 1603.6+2600 a.k.a. UW CrB, has a short P orb 
and likely shares properties with the other short period sys- 
tems, such as having a q <0.25, making tidal disc precession 
likely. 

3.3.8 Sco X-l (Fig\MJ\) 

The previously reported P SU p was incorrectly quoted as ~62 
d in OD01, instead of ~62 hrs (Kudryavtsev et al. 1989). 
Neither of these are significantly detected here. However, 
longer term periods ~200 d are strongly detected and appear 
to recur after essentially disappearing from MJD~51000- 
52000. 

Sco X-l is located in the OD01 intermediate instability 
zone on FigfT] where discs are expected to alternate between 
being warped and flat. However, the q ~0.7 used in OD01 
was only an estimate, and it has since been determined to be 
q ~0.3 (Steeghs & Casares 2002), making tidal disc preces- 
sion also likely. Recalculation using the OD01 method gives 

~1.5, meaning Sco X-l remains in the intermediate in- 
stability zone. 



3.3.9 XI 636-536 (Fio ISBj) 

The Psup ~46 d superorbital period of this quasi-persistent 
LMXB shows similar behaviour to that of KS 1731-260 in 
that it is strongly detected at MJD~53000-54000 when the 
source has made a transition to a lower state (Shih et al. 
2005). It does not appear to be stable during either state, but 
appears persistent and steadily evolving during the first half 
of the low state. The second part of the low state contains 
detections of shorter periods. 

The location of X1636-536 on Fig[T] makes formation 
of a warped precessing disc unlikely since it is below the 
r a dd = r curve. Warps are also expected to form during the 
persistent states of quasi-persistent X-ray binaries, but here 
the periodic signals are stronger during the low state. The 
quasi-periodic superorbital period detected in X1636-536 is 
believed to be a result of X-ray state changes, characterized 
by the anti-correlation between hard and soft X-rays (Shih 
et al. 2005). The optical and soft X-rays have been found 
to be correlated, as the latter is reprocessed in the accretion 
disc (Shih et al. 2011). 



3.3.10 X 1730-333 (Fig\EM 

This prototypical LMXB (it is the Rapid Burster in the 
globular cluster Liller 1 (Apparao & Chitre 1979)) shows an 
initially strong detection of the Psup ~217 d, which appears 
to remain steady up to MJD~51500. However, during that 
time there are a number of shorter term periods, that do not 
appear to simply be harmonics. A period of mostly ~100 d 
is strongly detected and steadily drifting until MJD~52500, 
when a marked evolution from shorter to longer periods be- 
gins. 

The superorbital periods in X1730-333 show evolution 
but are not steadily evolving (in the SMC X-l sense), neither 
are they unstable (in the Cyg X-2 sense). It rather appears 
to be a variety of steadily evolving periods, which initially 
appear to have some connection with each other. They may 
be indicative of the time-scale on which bursts recur, rather 
than periodicities. The location of X1730-333 on Fig[T] is 
unknown. 



3.3.11 XI 820-303 (Fig\£@ 

The Psup ~172 d in this LMXB is strongly detected and 
appears stable. Its detection is weaker during MJD~51200- 
52200 (but still significant) and simultaneously there are 
detections of shorter term periods during that time. The 
<-^86 d period is the first harmonic thereof, but the ~65 d 
period is not a harmonic. 

The location of X1820-303 on Fig[TJ makes it highly 
unlikely that any superorbital period in the system would be 
related to a precessing warped accretion disc. However, since 
it has been interpreted as a triple system (Chou & Grindlay, 
2001), the period would rather be connected with the M 
variations resulting from the effects of the third body and 
is therefore expected to be (extremely) stable. The minor 
variations visible in the DPS for the P 3 u P ~172 d are within 
conservative error estimates for the periods determined in 
each window, and are therefore not significant. 
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Table 3. Characterization of P SU p behaviour in HMXBs 



Table 4. Characterization of P S up behaviour in LMXBs 



Source 


Behaviour 


Mechanism 


Types^) 


Cen X-3 


varies 


mode 0-1* 


P,E,C 


Cyg X-l 


steady 


mode 0* 


U,R 


LMC X-3 


varies 


&M/r add = r * 




LMC X-4 


steady 


mode 0* 




SMC X-l 


varies 


mode 0-1* 




SS433 


steady 


jet/mode 0* 




X01 14+650 


steady 


BeX 


P,C? 


X1907+097 




mode 1+* 


P,T,C 


XTE J1716-389 


steady 


unknown 





* OD01 

M Liu, van Paradijs & van den Heuvel (2006) and ref. therein, 
P: X-ray pulsar, T: transient X-ray source, 
E: eclipsing system, R: radio emitting HMXBs, 
C: cyclotron resonance scattering feature at X-ray spectrum, 
U: ultra-soft X-ray spectrum 

3.3.12 XI 91 6-053 (Fig&ty 

There is an initial strong detection of the P S u P ~199 d before 
MJD~51000, but we found no evidence for any modulation 
near 5 d. Several periodic signals are detected initially, but 
none that remains stable or show steady evolution. 

As a double degenerate LMXB, the location of X1916- 
053 on Fig[T] makes it highly unlikely that any superorbital 
period in the system would be related to a precessing warped 
accretion disc. However, with q well below 0.25 (OD01) it 
would be very likely highly susceptible to disc precession 
due to tidal interactions with the donor (see e.g. Homer et 
al 2001). 

3.3.13 XI 957+ 115 fFig [2B} 

The results for this LMXB show detection of a P sup ~117 d 
initially and at MJD~52250, which appears to form part of 
an unsteadily evolving feature in the DPS. There are longer 
term features that are more significantly detected, but they 
are also unstable. 

According to the location of X1957+115 on Fig[T] it 
is not expected to develop a warped precessing disc. OD01 
suggest periodic signals may rather be a product of the time- 
scale of features in the lightcurve, due to changes in M. 

3.4 Summary of Results 

A summary of all 25 sources is presented in Table [3] & 2] 
for HMXBs and LMXBs respectively. Therein, concise de- 
scriptions of the behaviour of their superorbital periods are 
given, together with the mechanisms believed to be responsi- 
ble. Relevant references regarding mechanisms are contained 
in the text under discussion of each source. The tables also 
include the source types included in the catalogue papers by 
Liu, van Paradijs & van den Heuvel (2006,2007). 

Chaotic refers to unstable or highly variable periodic 
signals, like those displayed by Cyg X-2. Varying includes 
evolution of persistent periodic signals (like SMC X-l) or 
intermittent periodic signals (like X1636-536). Steady refers 



Source Behaviour Mechanism Types^ 



Cyg X-2 


chaotic 


mode 1+* 


B,Z,R 


EXO 0748-676 


varies 


A A I 


T,B,D,E 


GRS 1747-312 


varies 


unknown 


T,G,B,D,E 


GX 339-4 


varies 


AM 


T,U,M,R 


GX 354-0 


chaotic 


unknown 


B,A,R 


Her X-l 


steady 


mode 0* 


P,D,E 


IGR J17098-3628 


varies 


AM 


T,R? 


KS 1731-260 


varies 


AM 


T,B 


LMC X-2 






Z 


MS 1603.6+2600 






B 


Sco X-l 


chaotic 


intermediate* 


Z,M,R 


X1636-536 


varies 


AM 


B,A 


X1730-333 


varies 


unknown 


T 


X1820-303 


steady 


triple 


G,B,A,R 


X1916-053 


varies 


supcrhumps 


B,A,D 


X1957+115 


varies 


AM 


U 



* OD01 

t^l Liu, van Paradijs & van den Heuvel (2007) and ref. therein, 
A: known atoll source, B: X-ray burst source, 
D: dipping low-mass X-ray binary, 
E: eclipsing or partially eclipsing low-mass X-ray binary, 
G: globular-cluster X-ray source, 
M: microquasar, P: X-ray pulsar, 
R: radio loud X-ray binary, T: transient X-ray source, 
U: ultra-soft X-ray spectrum, Z: Z-type source 



to persistent or intermittent periodic signals that display 
no significant variation (like Her X-l & LMC X-4), i.e. the 
variations are within conservative error estimates for the 
periods. 



4 CONCLUSIONS 

This DPS overview aims to provide the long awaited time- 
resolved analysis that OD01 suggested would be required for 
aperiodic sources, as a basis for future investigations. An in- 
depth analysis of each source and detailed discussion of the 
full implications of these results, goes beyond the scope of 
this work. Here we wish to highlight whether the results from 
our DPS analysis suggest steady, evolving, persistent or un- 
stable periods, and how they relate to previously published 
periods. 

During the discussion of each source, we first compared 
the observed behaviour apparent from our analysis with 
the OD01 predictions for stability of accretion discs against 
radiation-driven warping, since the OD01 process has a de- 
pendence on system parameters (determining an XRB's lo- 
cation on Fig[T} which allows its feasibility for superorbital 
varying systems as a class to be assessed. Sources for which 
other mechanisms can account for the P Bup , were also com- 
pared to those predictions, where applicable. Such mech- 
anisms include tidal disc precession, precessing relativistic 
jets and M variations due to a third body, state changes or 
variations in the size of the equatorial disc surrounding a Be 
star donor. 
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4.1 Warped/Tilted Accretion Discs 

OD01 predicted that only a small fraction of X-ray binaries 
should display steady superorbital periods associated with 
stable, steadily precessing, radiation-driven warped accre- 
tion discs. Furthermore, they suggested that it would be 
even less common in LMXBs, since only those with P or t >1 
d are expected to produce warps. However, 16 of the 25 
sources contained herein are LMXBs. With 114 HMXBs and 
187 LMXBs contained in the catalogues of Liu, van Paradijs 
& van den Heuvel (2006, 2007), it implies that superorbital 
periods have been detected in a slightly larger fraction of 
LMXBs than HMXBs. 

However, recent Smooth Particle Hydrodynamics 
(SPH) simulations by Foulkes, Haswell & Murray (2010) 
on SMC X-l, Cyg X-l, Cyg X-2, X1916-053, LMC X-3, 
Her X-l, SS433 and a generalized LMXB, produced warps 
over all orbital periods considered. Contrary to OD01, they 
predict that superorbital periods should be very common 
in LMXBs, suggesting that the analytical OD01 approach 
is necessarily approximate and that the SPH simulations 
should incorporate the complexities involved in an irradi- 
ated accretion disc more accurately as their "disc continu- 
ously flexes in response to the changing orientation of the 
Roche potential". 

It is important to remember that the OD01 predictions 
for stability of accretion discs against radiation-driven warp- 
ing, are for specific a and e values and that the instability 
criteria depend largely on a and to a slightly lesser extent 
on e, making the predictions approximate. Different a and e 
values would yield very different predictions and these values 
are expected to differ from system to system. Other mech- 
anisms, such as wind-driven tilting and magnetic warping, 
may also produce warped/tilted accretions discs. 

Furthermore, warps need to be sustained for prolonged 
intervals to produce steady superorbital periods (as in Her 
X-l & LMC X-4). Lodato & Price (2010) used SPH simula- 
tions to consider the diffuse propagation of warps in viscous 
thin discs, where they determined the diffusion coefficient to 
be ~ 1/a (for small amplitude warps and a <0.1), in gen- 
eral finding that higher viscosity leads to slower diffusion 
and lower viscosity to more rapid diffusion. 

Of the 25 X-ray binaries considered here, 15 could be 
directly compared to their OD01 predictions. The results are 
generally consistent with their predictions. 

Firstly, our results support the OD01 prediction that 
LMXBs are unlikely to produce stable precessing warped 
discs, since they are mostly located below the r a dd = r line 
and their accretion disks are therefore too small to become 
unstable against warping in the expected radiation field. Her 
X-l remains the only LMXB to have a steady superorbital 
period that can be associated with its accretion disc, since 
X1820-303 is likely a triple system. Her X-l is not a typi- 
cal LMXB, since it has higher donor mass than most other 
LMXBs and therefore it might share some HMXB properties 
or be an Intermediate-Mass X-ray Binary (IMXB; see Pod- 
siadlowski, Rappaport & Han 2003). One of these properties 
appears to be its ability to produce a warped precessing disc, 
with its associated steady superorbital period. 

Secondly, OD01 predictions implied that HMXBs are 
more likely to produce stable warped precessing discs, al- 
though eccentric orbits would suppresses stable radiation- 



driven warping. Our DPS analysis has shown that not only 
LMC X-4 displays a stable superorbital period, but also Cyg 
X-l, SS433, X0114+650 & XTE J1716-389. Only the loca- 
tions of the latter two sources on Fig[T] are unknown. Cyg 
X-l, LMC X-4 & SS433 are located in the OD01 instability 
zone and therefore expected to have warped precessing ac- 
cretion discs that produce steady superorbital periods while 
they are persistent. However, the steady superorbital periods 
in SS433 and X01 14+650 may have other origins, since the 
former experiences relativistic jet precession and the latter 
is a BeX. 

Thirdly, we included 5 quasi-persistent sources for 
which superorbital periods have been published. Four of 
these are LMXBs (X1636-536, KS 1731-260, EXO 0748-676 
& IGR J17098-3628) and XTE J1716-389 is the only HMXB. 
As OD01 predicted, these quasi-persistent sources produce 
superorbital periods during their persistent (high) states, 
with the exception of X1636-536, where stronger detections 
occurred during the lower state. Its low state is certainly not 
an off state where the flux goes to zero, as seen in the other 
four sources. It also displayed an anti-correlation between 
hard and soft X-ray components (Shih et al. 2005). Further- 
more, only the HMXBs produced a relatively steady super- 
orbital period while the LMXBs all produced multiple peri- 
odic signals that were either evolving or unstable/chaotic. 

Finally, OD01 predicted that Sco X-l would be 
marginally unstable or stable and show variability, since it 
is located in the intermediate instability zone. Our results 
show long-term behaviour that may be interpreted in this 
manner. 

Clarkson et al. (2003) suggested that warps may man- 
ifest themselves as modulations in the X-ray flux, not only 
resulting from variations in the accretion disc structure but 
also by modulation of M. Varying absorption of X-rays from 
the central source by a warped inner disc or the variation 
in the uncovered X-ray emitting area apply to the former, 
while variation of the accretion rate onto the compact object 
or modulation of the mass flow rate through the disc (Dubus 
2003) apply to the latter. It was also suggested that compet- 
ing radiation-driven warping modes may cause variations in 
the warp itself, resulting in the steady evolution/ variation of 
the P 3U p in SMC X-l, consistent with its OD01 prediction. 
Her X-l, with P sup ~33-37 d (Leahy & Igna 2010) might be 
subject to the same effect, to a lesser extent. 

4.2 Complex Al variations 

Many of the superorbital periods already detected may be 
the result of the modulation of the M, rather than only be- 
ing associated with warped precessing accretion discs. Many 
may therefore also indicate the time-scale on which transi- 
tions occur between high and low flux states. The HMXBs 
Cen X-3, Cyg X-l & LMC X-3 and the LMXB X1957+115 
have very similar lightcurves that clearly show these tran- 
sitions on varying time-scales. However, in the case of Cyg 
X-l there is also evidence for a precessing warped accretion 
disc (Zdziarski, Pooley & Skinner 2011). 

In fact, the majority of the superorbital periods are 
clearly not steady, with many of those showing unsteady 
evolution of multiple periodic signals while they are persis- 
tent. Such behaviour may also be associated with a variety 
of time-scales produced by M variations. 
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4.3 The strength of time-dependent analysis 

Our results emphasize the importance of using a time- 
dependent periodic analysis approach such as DPS when 
investigating the long-term periods that have been as- 
sociated with accretion disc properties such as preces- 
sion/warping/tilting and M variations due to X-ray state 
changes. Traditional period analysis, which use the entire 
datasets at one's disposal, would miss a source like KS 1731- 
260 since its overall periodogram shows no significant peaks. 
Moreover, all the variations and complex behaviour of quasi- 
periodic signals go completely unnoticed. For example, the 
lightcurves of LMXBs X1730-333 & GRS 1747-312 show 
some similarity to Her X-l, but their superorbital periods 
are not steady and appear to be evolving or varying. 

Consequently it is not unexpected that studies of su- 
perorbital behaviour in hard X-rays (14-195 keV) using data 
from the Swift Burst Alert Telescope and soft X-rays (1.5-12 
keV) using RXTE ASM data, found no significant periodic 
modulations for X1916-053, Cyg X-2 & Sco X-l, while being 
able to confirm the previously reported periods in X1820- 
303 & X1636-536 (Farrell, Barret & Skinner 2009). Unstable 
and/or intermittent periodic signals may lie below the white 
and/or red noise level when considering the longer datasets, 
while the time-dependent DPS analysis shows intermittent 
and sustained periodic signals that appear to vary over time. 

Since we are no longer limited by desktop technology, 
applying this analysis method to the other XRB sources 
contained in the RXTE ASM is relatively straightforward. 
It may reveal complex periodic behaviour that has thus far 
gone unnoticed. 
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